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Abstract 

Restriction on the neutron-antineutron oscillation time in vacuum is obtained from 
latest SNO data on the deuteron stability, td > 3. 01' 10^^ years. Calculation performed 
within the quantum field theory based diagram technique reproduces satisfactorily results 
of the potential approach previously developed. The dependence of the obtained restric- 
tion on the total spin of the annihilating NN system and the deuteron wave function 
modifications is discussed. 
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1 Introduction 



Searches for the baryon number violating processes predicted by Grand Unified Theories (GUT) 
remain to be an actual experimental task during many years. The neutron- antineutron tran- 
sition induced by the baryon number violating interaction {AB — 2) predicted within some 
variants of GUT has been discussed in many papers since 1970 [1] , sec [2] — [7] . Experimental 
results of searches for such transition are available, in vacuum (reactor experiments [8], and 
references therein), in nucleus ^^O [9] and in Fe nucleus [10], in neutron magnetic trap [11], see 
[3, 12]. 

Here we derive a restriction on the neutron-antineutron transition time in vacuum using 
new data on the deuteron stability obtained in Sudbury Neutrino Observatory (SNO) [13] and 
our former result on the suppression of n — n transition in deuterium [14]. Restrictions on the 
n — n transition time in vacuum which follow from data on nuclei stability have been obtained 
previously in a number of papers by different methods, beginning with the potential approach 
of [4, 5, 6], see [15] - [21]. 

To introduce notations, let us consider first the nn transition in vacuum which is de- 
scribed by the baryon number violating interaction (see, e.g. [2, 16, 18]) V = ^nn(Ji/2, ai being 
the Pauli matrix, is the parameter which has the dimension of mass, to be predicted by 
grand unified theories and to be defined experimentally ^. As usually, a point-like n — n coupling 
is assumed here. The n — n state is described by the 2-component spinor ^, lower component 
being the starting neutron, the upper one - the appearing antineutron. The evolution equation 
is 

i'^^{Vo + V)^, (1) 

where the diagonal matrix Vq has matrix elements Vq^ — Vq^ — — i^nj"^ in the rest frame 
of the neutron (tbtv is the neutron (antineutron) mass, 7„ - the (anti) neutron normal weak 
interaction decay width, and we take ~ 7n, it follows from CPT-invariance of strong and 
weak interactions). Eq. (1) has solution 

^{t) = exp [-i {i^nnt ai/2 + Vot)] = 

Here is the starting wave function, e.g. for the neutron in the initial state = (0, 1)^. In 
this case we have for the wave functions at arbitrary time t 

^(n,t) = —isin ^^"' exp{—iVot), ^(n, t) = cos ^™'' exp{—iVot), (3) 

which describe oscillation n — n: In vacuum the neutron goes over into antineutron, also the 
discrete localized in space state, which can go over again to the neutron, so the oscillation 
neutron to antineutron and back takes place indeed. 

Since the parameter finn is small, the expansion of sin and cos can be made in Eq. (3) 
at not too large times. In this case the average (over the time t"^'* <^ 1/finh) change of the 

^ There is relation = '2,5m with the parameter 6m introduced in [2]. The neutron-antineutron oscillation 
time in vacuum is Tnn — I /5m = 2/ Unn, see also [21] and references in this paper. 



COS laism 



exp{-iVot)'i/o, 



(2) 



2 



probability of appearance of antineutron in vacuum is (for the sake of brevity we do not take 
into account the (anti)neutron natural instability which has obvious consequences) is 

,.2 j-obs 

W{n;t''''')/t''^' = \^{n,t''''')\yt''''' ~ ^^^^ (4) 

which has, obviously, dimension of the width F. From existing data obtained with free neutrons 
from reactor the oscillation time is greater than 0.8610^sec ~ 2.7 years [8], therefore, 

fXnn < 1.510-23 eV. (5) 



2 Suppression of n — n transitions in arbitrary nucleus 

Recalculation of the quantity fj^nn or Tnn from existing data on nuclei stability [9, 10, 13] is 
somewhat model dependent, and different authors obtained somewhat different results, within 
about one order of magnitude, see e.g. discussion in [18, 20, 21]. In the case of nuclei the n — n 
line with the transition amplitude ^nn is the element of any amplitude describing the nucleus 
decay A ^ (A — 2) + mesons, where {A — 2) denotes a nucleus or some system of baryons 
with baryonic number A — 2. The decay probability is therefore proportional to /x^^j, and we 
can write by dimension arguments 

FM ^ M -2) + mesons) ~ (6) 

mo 

where mo is some energy (mass) scale. 

It was argued in [14] that mo is of the order of normal hadronic or nuclear scale, mo ~ 
i^hadr ~ (10 — 100) MeV . The dimensionless suppression factor is therefore 

Fs = ^< 10-=^°, (7) 

mo 

We can obtain the same result from the above vacuum formula (4), if we take the observation 
time t°^^ ~ l/mhadr- 

The physical reason of such suppression is quite transparent and has been discussed in 
the literature long ago (see e.g. [2, 16, 19]): it is the localization of the neutron inside the 
nucleus, whereas no localization takes place in the vacuum case. 



3 The deuteron decay probability 

The case of the deuteron, which is most simple and instructive, can be treated using the 
standard diagram technique. Such technique or its modifications have been used in [17, 18] and 
[20] \ 

The point is that in this case there is no final state containing antineutron — it could be 
only the pn state, by the charge conservation. But this state is forbidden by energy conservation, 

^ Later the author of [17] tried to develope a field theory motivated approch to this problem which was 
criticized by scientific community [15, 16, 20, 14]. 
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since the deuteron mass is smaller than the sum of masses of the proton and antineutron. 
Therefore, if the n — n transition took place within the deuteron, the final state could be only 
some amount of mesons. The amplitude of the process is described by the diagrams of the type 




mesons 



Figure 1: The Feynman diagram describing n — n transition in the deuteron with subsequent 
annihilation of antineutron and proton to mesons. 



shown in Fig.l and is equal to 

, \ r. 2 f T(np mesons) d'^p , . 

T{D ^ mesons) = 2z^^„,m^/... j _ ^2 _ _ ^2 p 7^- (§) 



Here p and d are the 4-momenta of the virtual proton and deuteron. The constant gonp is 
normalized by the condition 



,2 



^ = — = ~ 0.049, (9) 

which follows, e.g. from the deuteron charge formfactor normalization = 0) = 1, see 

discussion and references in [14]. k = y/rriN^ — 46Met', en — 2.23 MeV being the binding 
energy of the deuteron. In the nonrelativistic reduction of Feynman diagram (Fig. 1) we should 
write for the vertex D ^ np 2mNgDnp and mNl^nn for the n ^ n transition amplitude, to ensure 
the correct dimension of the whole amplitude (here we correct some inaccuracy of our former 
consideration in[14]). 

Presence of the second order pole in the energy variable of intermediate nucleon is a 
characteristic feature of diagrams describing the n — n transition in nuclei, as discussed in [14]. 
This does not lead, however, to any dramatic consequences, and the integration over internal 
4-momentum d'^p in (8) can be made easily taking into account the nearest singularities in 
the energy po = in the nonrelativistic approximation for nucleons. The integral over d^p 
converges at small p k which corresponds to large distances between nucleons in the deuteron, 
r ~ 1/k. By this reason the annihilation amplitude can be taken out of the integration on the 
mass shell in some average point, and we obtain the approximate equality [14] 

T{D — )■ mesons) = 2gDnp^%IJ"nnlDNNT{np — )■ mesons) (10) 
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with the integral 



IDNN 



d^p 



(27r)4(2m)3 J - mN - P^/(2mAr) + i8] [mo - tun - Po - P^/(2mjv) - i5f 

d^p 1 



(27r)38mAr[fi:2+p^]2 647rmjv/«' ^^^"^ 

There is close connection between the amplitude of the n — n transition in deuteron and the 
deuteron charge formfactor at zero momentum transfer, which contains the integral Idnn as 
well [14], 

Using the standard technique, we obtain for the decay width (probability): 

V{D mesons) ~ l^inODnp^DNN''^^ j ~^ mesons)]^ d^{mesons), (12) 

^{mesons) is the final states phase space. Our final result for the width of the deuteron decay 
into mesons is 



D-^mesons 



-m 



N 



ann,S=l 



(np) 



rriN 



Pc.m.^np 



ann,S=l 



(13) 



where Pc.m. is the (anti)nucleon momentum in the center of mass system. This result is close to 
that obtained in [17] and in [18] where it was obtained using the induced np wave function. ^. 
The annihilation cross section of the antineutron with velocity vq on the proton at rest 

equals 



(T?p" = a{np — >■ mesons) 



1 



voa. 



np 



T{np — >■ mesons) I (mesons). (14) 
~ (50-55)m6 ~ (130-140) GeV'^ where 



According to PDG at small vq, roughly, 

averaging over antineutron and proton spin variables is assumed. We obtain from Eq. (13) for 
the deuteron life time 



Sttk 



IGttk 



/^nn mjv 



Pc.m.^np 



ann,S=l 



Pc.m.-^O 



' 2 

mif 



ann,S=l 
VoCTfip 



(15) 



Vq is the antiproton velocity in the laboratory frame where proton is at rest. 
If we define the suppression factor (dimensional) i?D as usually. 



(16) 



^The result of [18] given by Eq. (17) can be rewritten in our notations as 

^2 



D—ymesons 



0.01/x: 



nn 



K 



■uo-s-O ' 



(17') 



which is abot twice greater than our result. 
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with Tnn — 2//x„n, then we have from Eq.(15) that 



R 



IQttk 



D 



— ann,S=l 



(17)) 



Quite naturally suppression increases with increasing binding energy of the deuteron e^, or k. 
Numerically wc have Ru ~ 2.94'10^^sec~^ if we neglect possible dependence of the annihilation 
cross section on the total spin of the n — p system. This is in agreement with results of the 
papers [6] where Rd was found to be in the interval between 2.410^^ and 2.7510^^ sec~^ for 
different variants of the potential model. 

To get an idea what happens for heavier nuclei, we can use the formula (17) with greater 
value of e, about ~ 8 MeV, the binding energy per one nucleon in heavy nucleus. This leads to 
the value Ra ~ 5.5 10^^ sec^'^ sec^^, in qualitative agreement with [6]. We obtain same estimate 
if we replace the deuteron size ro = ^/^/^rriN by average internucleon distance in heavy nucleus, 

Tinter-N ~ 2 Fm. 

4 The role of the total spin dependence of the ftp anni- 
hilation cross section 



This result can be refined taking into account spin dependence of the NN annihilation ampli- 
tudes and going beyond the zero range approximation for the deuteron wave function. 

Within the deuteron neutron and proton are in triplet state, and after n — )■ n transition 
antineutron and proton remain in triplet state, therefore, only part of the total annihilation 
cross section corresponding to triplet state works in our case. 

We can write for the annihilation cross section averaged over spin states 



1 ann,S=l i 

where a^p and a^p 



-ann 
np 

ann,S=0 



ann,S=l , „ann,S=0\ 



np 



np 



3^ + 1 ann,S=l 



are the triplet and singlet annihilation cross sections, r — a.^ 



(18) 

,ann,S=\ i ann,S=0 



or cr, 



^ann,S=l 

np "np 

Prom these results and definitions we obtain the final formula 



np 



L / 



np 



4ra?"V(3r + l). 



rm 



'N 



td- 



np 



The supression factor is 



R 



D 



47r«;(3r + 1) 

4(3r + l)-nK 



1/2 



rm 



N 



np 



(19) 



(20) 



Numerically Rd ^ 2M-W^{2>r + l)/4r sec-^ and from new SNO data [13] td > 2>M IQ^^ 
years, we obtain 

Tnn > 1.8-10^ [4r/(3r + 1)]^^^ sec. (21) 
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At r = 1 we recover our former result for the spinless case, if r ^ 1 we would obtain 

Tnn> 2.1-10^ sec (22) 

In the paper [22] it was obtained from the combined analysis of the pp atom data and results 
of OBELIX scattering experiments at LEAR that for the antiproton-proton interactions rpp ~ 
0.42. The pp state is a mixture of isoscalar and isovector, whereas np is pure isovector, but if 
we take the same ratio for the np interaction, we obtain from the SNO data [13] that Tnn > 
1.55T0« sec. 



5 Possible deuteron wafe function modifications; 
conclusions 

To go beyond the zero range approximation for the deuteron wave function we can use e.g. the 
Hulthen deuteron wave function in the form 



1 ^aKja + k) [exp(-Acr) - exp(-Q;r)] 

^D,H[r - , [Z6 

V27r a - K r 

where a ~ 270 Mev is the Hulthen parameter [23]. In the limit o; — > oo we obtain the zero 
range deuteron wave function = ^K/2nexp{—K,r)/r. 

The additional formfactor appears in the integrand of the amplitude in Eq. (11) when 
we are using the Hulthen wave function. Instead of the integral 



d^p TT 



2 



[ac^+P^]^ k 
d^p 



^a{a + K)^ / r 2^ ^I2f 2^ = ~\/^~- (^4) 

^ J [K,^ + P^\^ [a^ + p^) ^ \l n 




we obtain now 



So, using more realistic deuteron wave function leads to additional factor — 0.92 which 

' ° y tt+re 

does not change substantionally our results: the final estimate for Tnn decreases by about 4%. 
This means also that possible modifications of the deuteron wave function at small internucleon 
distancies will not change our results considerably. 

To conclude, the results presented here arc in satisfactory agreement with previously 
obtained results of the potential approach by Dover, Gal and Richard [6], and also in crude 
agreement with [18, 17]. Our method opens a way to study possible relativistic effects and 
corrections. Restriction on the n — n transition time in vacuum obtained from the SNO data 
[13] is comparable with that from data on heavier nuclei stability. The best restriction which 
comes from data on oxygen stability, is Tnn > 3.3T0^sec [21]. 
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